Bound Modes in Dielectric Microcavities 
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We demonstrate how exactly bound cavity modes can be realized in dielectric structures other than 3d pho- 
tonic crystals. For a microcavity consisting of crossed anisotropic layers, we derive the cavity resonance fre- 
quencies, and spontaneous emission rates. For a dielectric structure with dissipative loss and central layer with 
gain, the 13 factor of direct spontaneous emission into a cavity mode and the laser threshold is calculated. 

PACS numbers: 42.60.-v,32.80.-t,78.67.-n. 
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I. INTRODUCTION 

One of the motivations for making a small and efficient 
laser is to enter a regime where quantum correlations in the 
emitted light can be created and studied. For this purpose 
spontaneous photons must be emitted into the cavity mode 
with high efficiency (the (3 factor), so that the noise from the 
random emission in other modes is small ||]. At the same 
time, the lifetime of the lasing mode needs to be large enough, 
to allow for stimulated emission. The fabrication of dielectric 
microstructures, like microspheres pL Hi and VCSELs (ver- 
tical cavity surface-emitting lasers)^^^, ^ ||] for the real- 
ization of strong coupling and high finesse is therefore a big 
challenge for cavity QED-type experiments. From a funda- 
mental point of view this also raises the question: what are 
the different types of systems where a small mode volume is 
still compatible with little leakage? Neither an ideal sphere 
[^, [l^] nor a stack of disks (an ideal VCSEL), for example, 
give rise to mathematically bound solutions. These systems 
have resonances with a finite width, which only vanishes for 
large dimensions. When these structures have dimensions of 
the order of one optical wavelength, the field strength of the 
mode may become large (and thereby the coupling constant), 
but at the cost of an increasing loss rate. Systems with ex- 
act bound states can exist in so-called photonic materials. The 
well known examples are photonic crystals with a point defect 
JTl] , and disordered structures [ [T3[ p4[ p^ . These systems 
are interesting for cavity QED, because ideally (no dissipation 
and infinite size) the lifetime of the bound state is infinite and 
propagating modes occur at other frequencies only. Photonic 
materials are complicated structures, however, so that fabrica- 
tion and modeling is not easy. In this paper we demonstrate 
the occurrence of exact 3d-bound states in dielectric structures 
that are not a photonic crystal. 



n. BOUND MODES IN ANISOTROPIC DIELECTRIC 



index of refraction. The occurrence of bound states relies on 
the use of dielectric material with particular polarization prop- 
erties. As we will show, the structure of Fig. |l] is in fact only 
one in a class of systems with bound modes. The systems in 
this class are characterized by the dielectric function eir) with 
the tensor form 



e{r)lei = 1 + zz{U{x) + y(y)] + (1 -zz)W{z). 



(1) 



Here e\ is the dielectric constant of the (isotropic) background 
medium. The functions U, V, and W represent the spatial de- 
pendence of the vertical and horizontal layered structures. In 
the case that the layers have a lower susceptibility than the 
background (for instance, when they represent air holes in di- 
electric material), U,V or W will be negative. These func- 
tions will be specified later. In a dielectric environment of the 
form (|l]), fields with a polarization in the horizontal plane will 
only be affected by the structures described by W{z), while 
vertically polarized fields are only perturbed by the structures 
described by U{x) + V{y). Moreover, because the structures 
depend on one or two spatial coordinates only, the propaga- 
tion of these particular fields in the other directions is free. 
Here we consider another specific type of polarization, which 
is governed by all three functions U , V , W . The electric field 
strength of this type can be expressed as 



Eir) = {k^y/k^ -zd/dz)f(x)g(y)h(z). 



(2) 



in terms of scalar functions fix), g(y), h(z) and eigenvalues 
k~, k. This is proven by direct substitution into the Maxwell 
equation V x V x £ = (uj/c)^e{r)E. After projecting out the 
Cartesian vector components, one finds that this equation sep- 
arates into independent equations for the scalar functions: 



-(d/dxffix) = 
-idldyfgiy) = 
-{dldzfhiz) = 



klf{x) + (k^-e)U{x)f{x), (3) 
l^.g(y) + {k^-k^^)V(y)g(y\ (4) 
k^Mz) + klW{z)h(z). (5) 



The system with the simplest design is shown in Figure ^ 
This cavity consists of horizontal layers reminiscent of a VC- 
SEL, with two crossed vertical planes of a material with lower 
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The separation constants kx, ky, and fc, in these equations, 
which denote wave vector components, must be related by 
k\ + k^. + k\ = k^ , with the length of the wave vector given by 
k = ej'^w/c. It can now be verified directly, that the field (^ is 
transversely polarized in the isotropic background. Note also 
that components normal to interfaces (where U,VoiW have 
steps) of the electric field E, the displacement field D = e(r)E 
and of the magnetic field B = -(i/uj)'V x E are all continuous. 
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That a dielectric structure of form with specific choice 
of U, V, and W can indeed support bound modes now follows 
from the following general considerations. Solutions that are 
localized in the x, y or z direction can be created when k^, or 
A:? is negative. Such solutions represent electromagnetic fields 
that are guided by vertical structures or a horizontal structure. 
Because Eqs. (^[^) are Schrodinger-type equations, it follows 
from standard wave mechanics that when U and V are posi- 
tive in a finite region, the prefactor - = k^. + kf, must be 
positive for such a guided wave. In that case both functions 
U{x) and V{y) act as attracting potentials, but only one of the 
two eigenvalues k^, k^ is negative. When, on the other hand, 
U and V both become negative, an attractive potential is ob- 
tained for negative k^ - k^, so that the sign in front of both 
functions is reversed. Therefore a wave can be localized both 
in the x and y directions simultaneously when U and V be- 
come less than -1/2. Eq. (||), is of the Helmholz type, with the 
effective potential W(z) multiplied by the eigenvalue. Only 
when W becomes less than -1, an attractive potential with a 
locaUzed solution occurs. Localized solutions for each of the 
scalar functions f,goih can also be obtained when the corre- 
sponding function U,V oiW is periodic in two half spaces. In 
that situation the wave function can be a defect state in a band 
gap. A fully bound state for the electric field will occur 
when all three scalar wave functions are localized simultane- 
ously. This can be achieved by means of combination of the 
different localizing effects. 

The conclusion of this analysis may be summarized as fol- 
lows. When the structures described by U and V have a higher 
index of refraction than the background, the functions U and 
V are positive and a 3d bound state requires periodic struc- 
tures in one of the two horizontal directions and also in the 
vertical direction, i.e. the system must be a 2d photonic crys- 
tal. When both U and V become less than -1/2, they describe 
crossed vertical slabs with a low index, and a 3d bound state is 
found with a Id periodic structure of horizontal layers. When 
W describes a horizontal layer of negative-index material, a 
Id photonic crystal of vertical layers is needed for a 3d bound 
state. Because k^ + ky + = k^, at least one of the eigenval- 
ues k^, k^, or k^ must be positive, so that localized solutions 
are not possible without the use of periodicity in at least one 
dimension. Because the type of systems proposed here is not 
a 3d photonic crystal, the frequency of a bound state will ly in 
the continuum of propagating modes. 

III. EXAMPLE: A ID PERIODIC STRUCTURE 

A. Mode Solutions 

We will now consider the most interesting situation: the 
Id periodic structure illustrated in Fig. [l|, which has a 3d 
bound state. The vertical and horizontal layers have widths 
d and b respectively and the layer-to-layer separation distance 
is a, so that the fundamental resonant wavelength lir/k^ will 
lie in the interval between 2a and 4-a. The anisotropy in- 
side the vertical layers is described with the dielectric tensor 
e = (l-zz)£i +zze2, while the tensor e = zzei + (1 -zz)e2 de- 




FIG. 1: Example of an empty-cavity geometry where bound field 
modes are realized. For the vertical and horizontal layers, an 
anisotropic dielectric structure of the form of Eq. (|l]) is required, as 
discussed in the text below. The layers are depicted transparent, be- 
cause their refractive index is lower than the index of the background. 
These layers, with their specific anisotropy, could be fabricated of air 
holes. The cylindrical boundary is not essential, because the bound 
modes decay exponentially with the distance from the origin. 

scribes the horizontal layers. The (real part of the) dielectric 
constant £2 must be smaller than e-i. This kind of anisotropy 
can be obtained by drilling air holes in isotropic high-index 
material in, respectively, the vertical and horizontal directions. 
We will discuss the cavity resonances of this structure and in- 
clude absorption in the layers by allowing an imaginary com- 
ponent of the dielectric function. The exact bound states of 
our system must then arize from the resonances in the limit 
of zero loss. This idea is illustrated in the thin-layer approx- 
imation, neglecting the overlap regions (intersections of two 
crossing layers) in order to obtain analytical expressions. 

The functions U and V in Eq. (|l]) describe the vertical low- 
index layers of width d and are expressed in terms of the pa- 
rameterx = {e2-S\)dle\. The horizontal layers with thickness 
b are described with W and are likewise expressed in the pa- 
rameter C = (£2 - S\)ble\. We adopt the forms 

U(x) = e{d-2\x\)xld, V{y) = e(d-2\y\)xld, (6) 
W{z) = XI 6'(/7-2|z-/fl|)e//7+6'(fo-2|z|)(a-0/^. (7) 

/=-oo 

where 0{x) the Heaviside function. The central horizontal 
layer at z = has a different parameter a = (£3 - e\)ble\. 
The situation a = describes the empty cavity of Fig. |l| We 
consider Im % > 0, Im ^ > to model dissipative structures 
and Im a < to represent a central layer with gain. It follows 
from continuity and the Bloch theorem that the desired solu- 
tions of Eqs. (^|-^) with potentials given by (^-^ are, outside 
the thin layers, of the form 

fix) = g(y) = e''>'^y\ 

hiz) = e'^'^^^P''smk,(\z\-la) + e"P''smk,(la + a-\z\), (8) 
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with I = int(|z|/a) the number of layers between position z and 
the origin. (Note that \z\-la and la + a- \z\ are the distances to 
the nearest layer below and above z.) The solution h{z) is an 
even function of z and for z > b/2 one has h(z + a) = e''"'h{z), 
so that p acts as a quasi momentum. Because the system is 
open and x, ce are not real, kx, ky, k, and p will generally 
be complex valued. The condition for an outgoing resonance 
implies Re k^ > 0, Re ky > 0, and Re p > 0. These wave 
vector components can be expressed in terms of k, and the 
material constants by 

kx = ky = 2i/(-2x-d), 

k^ = kl - 8/(2x + df. (9) 

e'P" = cos k,a - (k,al2) sin k,a. 

The values of k, are the solutions of the closed equation 

i2lk,f[\ + ^k, cotk.a] = ^2 _ _ ^)2^ (^Q-) 

Modes that decay in time (Im k < 0), describe loss of energy 
both by leakage of light out of the system to infinity and by 
absorption of the light in the structure. The amplified modes 
(Im A; > 0) are not realistic for long times, because saturation 
effects can not be neglected at high intensities. The spatially 
localized solutions describe the situation that the energy cre- 
ated at z = is fully dissipated in the other layers. In the limit 
of a passive cavity without loss, one has a = Im x = Im ^ = 0, 
so that the right-hand side of Eq. ( |l0| ) disappears. Then the 
k^ are real and ly inside the energy gaps of the band structure 
for the fully periodic potential W(z) (the case a = ^). Be- 
cause kx and ky are purely imaginary and k~ is real, k is real. 
This identifies the 3d-bound states. Outside the layers the field 
is exponentially decaying in all directions. The decay in the 
z direction in Eq. (|^) derives from the quasi momentum of 
the form p„ = mria + iq„ (the index n = 1,2, .. . labels the 
resonances). The solutions to lowest order in ^ are given by 
k^„ = rniia - and q„ = {mi^)^l2a^. 

We give two numerical examples for the geometry d = a = 
4b. For £i = 13, £2=1 (air holes in GaAs), one finds a period 
of a = 0.46A, (A is the wavelength of the first resonance) and 
qia = 0.29. When e\ = 2.3, £2 = -20 (silver in glass), then 
a = 0.68A, ^ia = 0.80. 



B. Spontaneous emission factor and laser threshold 

In presence of absorption, the imaginary part of the fre- 
quency uj„ = ckn/e\'^ of an outgoing resonance determines 
the decay rate of the probability that a photon is present in 
the corresponding mode. For small values of the absorption 
in the layers and for small gain, the following approximation 
applies. Define from here on the real and (small) imaginary 
parts of the effective 2d susceptibilities and the wave vector 
with the notation x + "7> £, + 'C and k„ - ij„ . The small decay 
rate of a resonance is then cyjel'^. This can be interpreted as 
an imaginary correction term to the real frequency ut,, of the 
3d-bound state. The result of the expansion of Eqs. (^-|l^ is 
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k„V2 {£, + a + k%£,^a)k„' 



(11) 



FIG. 2: Spontaneous emission rates I^, and Tp^op into a cavity reso- 
nance and into the propagating modes, for a horizontal dipole emitter 
in the center of the geometry of Fig. as a function of the wave vec- 
tor (the frequency is lj = ck/e\'-). The adopted effective 2d dielectric 
susceptibilities are x = -12a/13, ^ = ~3a/13, with absorptive parts 
77 = O.Olx, C = 0.01^. This could correspond lo d = a = 4b and a 
background dielectric constant ei = 13 and with £2 = 1- The dashed 
and dotted curves are the emission rates in the back ground and in 
vacuum. The width of the sharp resonances (n = 1,2,3) vanishes 
for = and when the absorptive loss is compensated by gain in an 
active layer at z = 0. 



Before including a linear gain medium at z = 0, we first con- 
sider spontaneous emission rates from a single dipole placed 
at the origin of the cavity. The probability that a sponta- 
neous photon is emitted directly into a cavity resonance is 
/3„ = r„/(r„ + Fpiop) in terms of the partial rate into a cavity res- 
onance r;, and into the total of propagating modes Fp,op. The 
contribution of the resonances to the emission rate is given by 
a sum of the Lorentzians 

r„(k) = (6nr/R„)j„/[(k-k„f + -fll 

Here, Rn = {1/^ + kl„a/ sin^ k,„a)(4x^ -d^)(k„/k,,y is the in- 
tensity of the (normalized) bound state at the origin, F = 
I2k^ l.i^ef^/Pi{2 + S])^ is the emission rate in the background 
medium [|l6|, and the width 7,, is given by Eq. ( [TTI ) for a = 0. 
The rates are plotted in Fig. ^ as a function of the wave 
vector for a fixed value of the effective susceptibilities x^ 
The efficiency factor at the resonances is plotted in Fig. 
||a). There are several types of propagating modes; their clas- 
sification and specific features will be discussed elsewhere 
Ip^]. For strong coupling, cj„/s\'^ « F, the decay is no 
longer exponential [ p^ and one enters the cavity QED regime. 
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FIG. 3: a) Spontaneous emission factor (3 into the cavity mode versus 
the dielectric constant of the structure £2, at the first three resonant 
values k„ (n = 1, 2, 3). These correspond to the maxima of the peaks 
in Fig. ^. b) Gain constant a versus dielectric constant ei that is 
required to reach the lasing threshold, according to Eq. (p^. The 
absorptive components of the effective dielectric susceptibility were 
taken as r; = O.Olx, C = 0.01^. 



For the parameter Im a > our model represents a cavity 
with a gain layer at z = 0. This gain compensates the loss in 
the layers and reduces the width of the cavity resonances, as 
can be seen from Eq. (11 ). When the loss is fully compen- 

arizes. This corresponds 



sated, the stationary situation 7,, 
to the lasing threshold. The gain a needed to compensate for 
a small loss, is linear in the loss constants 77, From Eq. ( pi] ) 
it follows that the required value of a is given, including the 
effect of the beta factor, by 



a = 



I3„ 



2i r C K-I&'^V ^ + a + ki„ea 



V2 



(12) 



The behavior as a function of the susceptibility ^ is plotted in 
Fig. ^p). We conclude that the required value of the effective 



gain constant a can easily be reached for semiconductor ma- 
terial (for which the value a = -IO/77 is quite realistic). The 
behavior of the laser above threshold is not discussed here, 
because then the stationary state is determined by saturation, 
not included in our model. 



IV. CONCLUSION 

Apart from absorption, a structure of finite size (layers of fi- 
nite extent and finite in number) will also give rise to leakage. 
This, of course, is also the case for a defect state in a common 
three-dimensional photonic crystal. The finite size approxi- 
mation is likely to be very good, however, since the (analytic) 
mode functions are exponentially localized evanescent waves. 
Therefore, the coupling to the outside world is exponentially 
small in the system size, while the mode volume remains of 
the order of one cubic wavelength. For a cubic system of size 
Nd X Nd X Na, the fraction of mode volume outside the system 
is roughly exp(-2A^|A:j|t/-A^^ifl). For > 6 (GaAs case), and 
for > 23 (glass/silver case), this fraction is less than 10"'". 
In this respect, the proposed device of Figure |l] provides an 
interesting alternative to the use of whispering gallery modes 
in dielectric spheres, where smaller mode volume leads to in- 
creasing loss. Semiconductor structures with 3d-bound states 
may be very promising for future cavity QED experiments and 
quantum communication applications. 
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